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ABSTRACT 

We develop a transfer matrix formalism for two-dimensional pure gravity. By 
taking the continuum limit, we obtain a "Hamiltonian formalism" in which the 
geodesic distance plays the role of time. Applying this formalism, we obtain a uni- 
versal function which describes the fractal structures of two dimensional quantum 
gravity in the continuum limit. 



Recent developments of two-dimensional gravity have provided us with an un- 
ambiguous definition of quantum gravity. This is based on the equivalence of the 
continuum formulation [1] and the dynamical triangulation. [2] In two dimensions, 
we thus have a regularized quantum gravity which has a definite continuum limit. 
The remarkable success of the matrix models [3] further elucidated topological as- 
pects of the theory. However we still lack a general formulation for describing 
quantum fluctuations of space-time and for evaluating physical observables such as 
fractal dimensions. In this paper we propose a new formulation which is a kind of 
Hamiltonian formalism for quantum gravity. We show that a geodesic distance de- 
fined on a dynamically triangulated surface can be regarded as the "time" variable 
for defining the transfer matrix. We then obtain a "Hamiltonian" in the continuum 
limit, and analyze the fractal structures of the space-time. 

Let us consider a cylinder with an entrance loop(c) and an exit loop(c'). (See 
Fig. 1.) We introduce the following quantity which is formally defined in the 
continuum framework: 



N(L,L';D;A) = f —^l-—^jd 2 x^-A)5{j y/g^dxvdx" - L) 

c 

8( ( yj g^dxv-dx" - L') J] 5(d(P, c)-D). 

(1) 

Here in the path integral the total area of the surface and the lengths of c and c' 
are constrained to A, L and L', respectively. We further impose a constraint that 
any point on d has a geodesic distance D from the entrance loop c. To be precise, 
the geodesic distance d(P, c) is defined as the minimum distance between the point 
P and a point on the loop c. Note that the exit loop c' is one of the loops which are 
composed of points having geodesic distance D from the entrance loop c. For later 
convenience we introduce a making point on the exit loop as is shown in Fig.l. 
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An important property of N(L, L';D;A) is the following composition law: 

oo oo 

N(L, L'\ D;A) — J dA' J dL"N(L, L"\ D'\ A')N(L", L'\ D-D 1 - A- A'). (2) 
o o 

In other words, if we define N(L, L'\ D) by N(L, L'\ D) = J™ dAN(L, L'\ D; A)e~ tA J 
where t is the cosmological constant, it can be regarded as the matrix element 
< L | e~ DH | L' > for a "Hamiltonian" H. In this sense we call N(L, L'\ D) 
the proper time evolution kernel. Roughly speaking, eq. (2) asserts that a cylinder 
with height D can be decomposed into two cylinders with height D' and D — D' . 
More strictly, there could be several loops which have geodesic distance D' from 
the entrance loop c. From only one of them, however, the exit loop c' has geodesic 
distance D - D' . (See Fig.2) 

We now give a constructive definition of N(L, L'; D; A) in terms of dynamical 
triangulation. The strategy to define its lattice counterpart N(l, I'; d; n) is the 
following. We first introduce the notion of deforming a loop one step forward, 
whose precise definition is given in the next paragraph. As we will see there, a 
loop can split into several loops after a one-step deformation. In that case each 
loop is independently deformed in the next step. By repeating this procedure d 
times, we have a rf-step deformation of a loop. Then N(l, I'; d; n) is defined as the 
number of possible triangulations of a cylinder satisfying the following conditions 
which correspond to the four (^-functions in eq.(l): (i) The number of triangles is 
n. (ii) The entrance loop (c) consists of / links, (iii) The exit loop (c ; ) consits of 
V links, (iv) The exit loop is one of the connected components obtained after a d 
step deformation of the entrance loop. 

In order to give a precise definition of the deformation, let us consider a loop c 
on a triangulated surface. To deform c one step forward means to remove triangles 
attached to c in the forward direction and two-fold links on c. A typical example 
is illustrated in Fig. 3, where c is represented by a solid line, and the forward 
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direction is assumed to be inward. In this example, c has three two-fold links on it 
as indicated by a, (3 and 7, but we still regard it as a single loop. After a one-step 
deformation, this c will split into three loops which are indicated by the dotted 
lines. The notion of two-fold links might seem an idle complexity. As we will 
see, however, it plays an important role in the explicit evaluation of the transfer 
matrix. Note that the deformation of a loop considered here is closely related to 
the geodesic distance on the dual lattice. Actually, if a loop d is obtained from 
loop c after a d-step deformation, any point on c' has geodesic distance d from c. 
This is why we impose the condition (iv) in order to express the last 5-function in 
eq.(l). 

Since N(l,l';d;n) is the lattice counterpart of N(L, L'; D; A), it satisfies the 
same composition law as eq.(2). Therefore, if we define the d-step evolution kernel 
N(l, I'; d) by N(l, I'; d) = J2n=o ^(1, 1'] d; n)K n and regard it as the /, V element 
of a matrix N(d), it can be decomposed into a product of single step evolution 
kernels: 

N(l,l';d) = (N(d)\ t , = (iV(l)V, (3) 

which means that N(l) plays the role of a transfer matrix. 

Next we show that the transfer matrix N(l) can be evaluated by a combina- 
torial analysis. To proceed, we first need to evaluate the disk amplitude 

00 

F(y;K) = ^^F(/;n), (4) 

l,n=0 

where F{1; n) is the number of possible triangulations of a disk which has a bound- 
ary of length I and consists of n triangles. The disk amplitude F(y; K) can be 
evaluated by the large- iV 3 matrix model [4] as 
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where a, b and c are functions of K determined by 



a=±-c+^2c(±-c), 

^ = c( ^- c)(c -2^ } - 
( b < < a < c ) 



(6) 



As is easily seen from (5) and (6), F(y;K) is analytic in both K and y with 
finite convergence radii. At the critical point where F(y; K) becomes singular as 
a function of K, the values of a and c coincide and the convergence radius with 
respect to y is equal to 1/a. These critical values can be easily obtained from eq.(6): 
K 2 C = 1/12^3, y c = (3V4 -3~V*)/2, a c = c c = l/y c , b c = (l-y/3)/2K c . 

Near this critical point we take the continuum limit by parametrizing K and 

y as 

K = K c e- eH , y = y c e~< . (7) 

Then a, b and c are calculated from eq.(6) up to order e as 



4 r 2 r 

a = a c - ^T/I e ^' b = & c, c = c c + ^YTjeV t. (8) 



We point out that the parametrization of eq.(7) is a natural one in the following 
sense. If we replace K n and y l in eq.(4) with K n = K™e~~ ne 1 and y l = y l c e~ leC > , 
respectively, it is clear that the continuum limit e — > corresponds to taking the 
n — > oo and / — > oo limit with the physical area A = ne 2 and the physical length 
L = le kept finite. Thus t and £ are the conjugate variables to the area and the 
boundary length, respectively. 
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Substituting eqs.(7) into eq.(5), we obtain the following disk amplitude near 
the continuum limit: 

P iC , t ) = i±^(l - V3<) + <i±^/(C, r)S + 0(0, (8) 

where 

/(C,r) = (2(-v^)\/( + V^, ^ = YT^ 1 ' (10) 

Here f(C,T~) is a universal function in the sense that it does not depend on the 
details of the triangulation prescription. 

We now use combinatorics to evaluate the generating function of the matrix 
element of the transfer matrix N(l): 

oo oo 

N(y,y';K) = 5>V''(#(1)) V = E vV K n N{l,l' -l-n). (n) 

1,1' l,l,n=0 

What we are going to do is to sum up all possible triangulations of a cylinder such 
that the exit loop d is one of the connected components of the one-step deformation 
of the entrance loop c. In order to define a natural matrix multiplication of the 
transfer matrix, we introduce a marking point on d but not on c.(See Fig. 4.) Let 
us first consider the triangle which is attached to the marked link. There are only 
three possible types for this marked triangle; 1), 2), and 3) in Fig.5a, where F in 
3) denotes an arbitrary triangulation of a disk which is connected to the marked 
triangle through a vertex. After a one-step deformation, c proceeds to the marked 
link on d and the disk part will be disconnected. As we can see in eq.(ll), the 
power of y, y', and K counts the number of links on c, links on d, and triangles 
respectively. Then the contribution to N(y, y'; K) from the marked triangles is 
easily evaluated as follows: 1) yy' 2 K, 2) yy' 2 K, 3) y 2 y'KF(y; K). 
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We then examine the structure of the rest of the triangulation by going along 
the exit loop c' . There are four different types of basic structures 1),2),3) and 
4) as depicted in Fig.5b. Their weights are easily evaluated as 1) yy' 2 K, 2) 
y 2 y'KF{y; K), 3) y 2 F{y; K), and 4) yK{F{y; K)-l). Starting with the marked 
triangle, we can attach any of these four structures one by one repeatedly, and then 
come back to the original marked triangle. Therefore the contribution of these four 
structures to N(y,y';K) can be expressed by a geometric series. We thus obtain 
the following form for the generating function of the transfer matrix elements: 

N(y,y';K) = {2yy' 2 K + y 2 y'KF(y; K)} 

oo 

x J2iyy' 2K + y 2 y' K F(y, k) + y 2 F{ y - k) + y K{P{ y] k) - i)} n 

n=0 

. 2yy' 2 K + y 2 y'KF(y;K) 

1 - yy' 2 K - y 2 y'KF{y- K) - y 2 F(y; K) - yK(F(y; K) - 1) ' 

(12) 

We next consider the continuum limit of this transfer matrix. In order to see 
how the continuum limit should be taken, let us express the composition law in 
terms of the generating functions: 

^j^N(y,z-dr,K)N{- z ,y'-d2-,K) = N(y, y'; d\ + K), (13) 

where N(y,y';d;K) = E/"o E"=o vW K n N{l, d;n). From eq.(13) it is 

clear that the continuum limit for y and y' should be taken around different values 
which are inverse to each other. Furthermore, since the structure of the entrance 
loop is similar to the boundary loop of the disk amplitude, it is natural to expect 
that the continuum limit can be taken by setting 

y = y c e-<, y' = y^e^' , K = K c e^ h . (14) 
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If this is true, the composition law (13) can be re-expressed as 



e 
2rri 



zoo 

: J dZN((,Z;d l] t)N(-Z,C';d 2 ;t) = N((X';di + d2;t), (15) 



where N(C, C'i d; t) stands for N(y, y'\ d; K) for the values of y, y' and K given by 
(14). 

The validity of this continuum limit is explicitly checked by substituting (14) 
into eq.(12). We obtain 

N(y,y';K) = N(C,?;t) = 1 ; + O(e ), (i 6 ) 

e C + C - ae2/(C, r ) 

where a = ^2/(9^-1) and f((, r) is given by eq.(10). In this eqation we have 
the following miraculous cancellations, which convince us that the continuum limit 
we consider here is in fact the right one. First of all 0(e°) terms in the denominator 
of N(y,y'; K) cancel out. Secondly, the coefficients of ( and (' are 1. Thirdly, the 
residue of N(y, y'; K) with respect to the C + C' term is 1/e, which, as we see below, 
exactly cancels the factor e in (15). These are highly nontrivial results. 

By substituting (16) into (15) for d\ = 1 and d 2 = d, we obtain 



N((,(';d+l;t) = -L / df- 1 r— - N(~t, d; t) , 



— ICO 



C + C- «e2/(C,r) (17) 
= N(C-aet/(C,r),C';d;0> 

which leads to 

7V(C,C';d+l;*) - N((,(';d;t) = - ay /tf(C,T)^N(C,?;d;t). (18) 

We then take the continuum limit of eq.(18) by introducing D = oi\fed and taking 
e — > 0. We thus obtain the following continuum differential equation for D) = 
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N((,(';D;t): 



^(C,D) = -f((,T)^((,D). (19) 



In other words, the proper time evolution kernel can be identified with the following 
matrix element: N((, C'; D; t) =< ( \ e~ DH \ (' >, where if is a "Hamiltonian" 
given by 

H = /(C,r)^. (20) 

It should be noted that the miraculous cancelations that occurred in eq.(16) are 
crucial in the derivation of the continuum differential equation (19). 

The initial value problem for the differential equation (19) can be easily solved 
by the use of the characteristic curve method as ^((;D) = ^((';0), where (' is 
defined by f^, d("/ f((", r) = D as a function of ( and D. We can now evaluate 
the proper time evolution kernel N(L, L'; D) which is given by the inverse Laplace 
transformation of N(C, C'j D; r) for ( and By carrying out an inverse Laplace 
transformation for the solution of (19) with the initial condition tp(C', 0) = e~ L *>, 
we obtain 



•too 

N(L,L';D) = ^- J dCe^e" 1 ^'. (21) 



Note that eq.(21) indeed reproduces the corresponding initial condition N(L, L'; 0) =| 
S(L — I/), because we have (' = ( for D = 0. 

For later use we give N(L, L'\ D) for small values of L and r: 



L' r 2 3 ,„ L' , 

/ 7fZ' 

+ 0(L°). 

(22) 



* N(L, L'; D) is related to its discrete version N(l, I'; d) by N(L, L'; D) = e^V N(l, I'; d). 
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We also calculate the inverse Laplace transformation of the disk amplitude as a 
double series expansion with respect to e and t: 



too 

F(L) = ±, J d( e L < F(C, f)e- 3/2 



= ^(i-Viie)^ (23) 

+ ! ^ (l ^_^ + ^ )+0(A 

where F((;t) is given by eq.(9). 

We now apply the formalism developed here to the analysis of the fractal 
structures of the space-time. [5, 6, 7] We consider a large enough space-time with 
spherical topology for two-dimensional pure gravity. We take an arbitrary point 
P and consider the set of points S(P; D) whose geodesic distances from P are less 
than or equal to D. The boundary of S(P; D) usually consists of many loops with 
various lengths. Let p(L; D)dL be the number of loops belonging to the boundary 
of S(P; D) whose lengths lie between L and L + dL. As we explain below, the 
quantity p(L; D) can be evaluated as 



lim &mLoL.,D )iF(L) 
^O'"" &i,F(L«) 



(24) 



where 



G(x) = Y^(*- 5/2 + \^ 12 + y* 1/2 )e"*, (25) 



with x — L/D ' as a scaling parameter. In eq.(24), the limit Lq — > corresponds 
to shrinking the entrance loop to a point P, and the limit r — > is taken to have 



10 



the thermodynamic ilmit. The second derivatives with respect to r in eq.(24) are 
introduced to avoid small area dominance in the r — > limit. Any higher derivative 
works for this purpose without changing the final answer. The numerator on the 
RHS corresponds to an operation of gluing a disk with a boundary of length L 
to the exit loop of the cylinder. The factors 1/Lq and l/L in front of F's are 
introduced to convert marked boundaries to unmarked ones. As we have already 
explained, there are several loops which have geodesic distance D from the entrance 
loop. Recognizing that the exit loop is one of these loops, we easily realize that 
the ratio in eq.(24) counts the number of loops which have the boundary length L. 

One of the surprising properties of the function p{L\ D) is that it is essentially 
a universal scaling function of the scaling parameter x = L/D 2 . For small values 
of L, however, p(L; D) includes non-universal parts which have negative power 
dependence on the lattice constant e. In order to examine the scaling property, it 
is convenient to introduce the following quantities: 

oo 

<L n > = J dL L n p(L; D). (26) 
o 

From eq.(24) it is easy to show that 

< L° > oc const x D 3 e~ 3 / 2 + const x De~ 1 / 2 + const x D°, 

< L 1 > oc const x dV 1 / 2 + const x D 2 , (27) 

< L n > oc const x D 2n (n > 2). 

Since D 2 p(L; D) is essentially the universal function G(x) of the dimensionless 
parameter x, it is an ideal quantity to measure in computer simulations. In fact a 
numerical study for pure gravity shows a clear universal scaling behavior according 
to G(x).[8] As we can see in eq.(27), < L° > and < L 1 > include e dependent 
non-universal part as a dominant contribution while < L n > for n > 2 includes 
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only universal part and thus is expected to show clear fractal behaviors. The con- 
troversial numerical results for the fractal dimension by Agishtein and Migdal [6] 
can be understood by the non-universal behavior of the corresponding quantities. 

Although we have only treated pure gravity (c = 0) in this paper, it is im- 
portant to investigate quantum gravity coupled to matter fields. In the limit 
c — > — oo it is known that the two-dimensional gravity system approaches to 
the classical limit. We thus expect that the wildly branching space-time sur- 
face approaches to a smooth surface. On the smooth surface there is only a 
single boundary at a geodesic distance D from a given point. We then have 
p(L;D) =5(L — 2-kD) = 1/DS(L/D — 2tc), which again suggests a scaling function 
with respect to a scaling parameter x = L/D. It is thus very natural to expect that 
the gravity, for example, with c = — 2 matter shows the similar scaling behavior. 
In fact our recent numerical investigation shows a clear scaling behavior of the 
function p(L; D) which has a similar scaling parameter x = L/D 19±0 - 2 but has a 
slightly different power behavior from (25). [9] It would be very interesting if we 
can find an analytic formulation for deriving p(L; D) in a model with c 7^ 0. 

In this paper we have given an explicit form of "Hamiltonian" by eq.(20) which 
could be viewed as a time evolution generator of closed string, where time is iden- 
tified with the geodesic distance. This formulation may provide a new formulation 
of closed string field theory. 
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FIGURE CAPTIONS 



Fig. 1 Cylinder with entrance loop c and exit loop d . 

Fig. 2 Decomposition of a cylinder with height D into two cylinders with height 
D' and D - D'. 

Fig. 3 A Typical example of the deformation of a loop. The dotted loops are 
obtained after a one-step deformation of the solid loop. 

Fig. 4 A triangulation contributing to N(y,y'; K). The dotted loop (d) is one of 
the components of the one-step deformation of the solid line (c). 

Fig. 5 a) Three types for the triangle attached to the marked link. b)Four basic 
structures which may appear along c. The solid and dotted lines stand for a 
part of the entrance(c) and the exit(c') loops, respectively. 
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Fig.3 




Fig.4 



